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Abstract 

In  proposals  for  quantum  computers  using  arrays  of  trapped  ultracold  polar  molecules  as  qubits,  a  strong 
external  field  with  appreciable  gradient  is  imposed  in  order  to  prevent  quenching  of  the  dipole  moments 
by  rotation  and  to  distinguish  among  the  qubit  sites.  That  field  induces  the  molecular  dipoles  to  undergo 
pendular  oscillations,  which  markedly  affect  the  qubit  states  and  the  dipole-dipole  interaction.  We  evaluate 
entanglement  of  the  pendular  qubit  states  for  two  linear  dipoles,  characterized  by  pairwise  concurrence,  as 
a  function  of  the  molecular  dipole  moment  and  rotational  constant,  strengths  of  the  external  field  and  the 
dipole-dipole  coupling,  and  ambient  temperature.  We  also  evaluate  a  key  frequency  shift,  Am,  produced  by 
the  dipole-dipole  interaction.  Under  conditions  envisioned  for  the  proposed  quantum  computers,  both  the 
concurrence  and  Am  become  very  small  for  the  ground  eigenstate.  In  principle,  such  weak  entanglement 
can  be  sufficient  for  operation  of  logic  gates,  provided  the  resolution  is  high  enough  to  detect  the  Am  shift 
unambiguously.  In  practice,  however,  for  many  candidate  polar  molecules  it  appears  a  challenging  task  to 
attain  adequate  resolution.  Simple  approximate  formulas  fitted  to  our  numerical  results  are  provided  from 
which  the  concurrence  and  Am  shift  can  be  obtained  in  terms  of  unitless  reduced  variables. 


*  Corresponding  email:  dherschbach@yahoo.com 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

28  FEB  2011 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2011  to  00-00-2011 


4.  TITLE  AND  SUBTITLE 

Entanglement  of  polar  molecules  in  pendular  states 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


6.  AUTHOR(S) 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROIECT  NUMBER 


5e.  TASK  NUMBER 
5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Department  of  Physics, Texas  A&M  University, , College  Station, TX, 77843  REPORT  number 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

In  proposals  for  quantum  computers  using  arrays  of  trapped  ultracold  polar  molecules  as  qubits,  a  strong 
external  field  with  appreciable  gradient  is  imposed  in  order  to  prevent  quenching  of  the  dipole  moments  by 
rotation  and  to  distinguish  among  the  qubit  sites.  That  field  induces  the  molecular  dipoles  to  undergo 
pendular  oscillations,  which  markedly  affect  the  qubit  states  and  the  dipole-dipole  interaction.  We  evaluate 
entanglement  of  the  pendular  qubit  states  for  two  linear  dipoles,  characterized  by  pairwise  concurrence,  as 
a  function  of  the  molecular  dipole  moment  and  rotational  constant,  strengths  of  the  external  field  and  the 
dipole-dipole  coupling,  and  ambient  temperature.  We  also  evaluate  a  key  frequency  shift,  4!,  produced  by 
the  dipole-dipole  interaction.  Under  conditions  envisioned  for  the  proposed  quantum  computers,  both  the 
concurrence  and  4!  become  very  small  for  the  ground  eigenstate.  In  principle,  such  weak  entanglement  can 
be  sufficient  for  operation  of  logic  gates,  provided  the  resolution  is  high  enough  to  detect  the  4!  shift 
unambiguously.  In  practice,  however,  for  many  candidate  polar  molecules  it  appears  a  challenging  task  to 
attain  adequate  resolution.  Simple  approximate  formulas  fitted  to  our  numerical  results  are  provided  from 
which  the  concurrence  and  4!  shift  can  be  obtained  in  terms  of  unitless  reduced  variables. 


15.  SUBIECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

Public  Release 

40 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


I.  INTRODUCTION 


Since  the  original  proposal  by  DeMille  [1J],  arrays  of  ultracold  (<  1  mK)  polar  molecules  have 
come  to  be  considered  among  the  most  promising  platforms  to  implement  a  quantum  computer  [2- 
13fl .  His  proposal  describes  a  complete  scheme  for  quantum  computing  using  as  qubits  the  dipole 
moments  of  diatomic  molecules,  trapped  in  a  one-dimensional  optical  lattice,  partially  oriented 
in  an  external  electric  field,  and  coupled  by  the  dipole-dipole  interaction.  The  qubit  states  are 
individually  addressable  because  the  field  has  an  appreciable  gradient  so  the  Stark  effect  is  different 
for  each  location  in  the  array. 

A  subsequent  proposal  has  advocated  coupling  polar  molecules  into  a  quantum  circuit  using 


superconducting  wires  [14] .  Such  capacitive,  electrodynamic  coupling  to  transmission  line  res¬ 


onators  is  analogous  to  coupling  to  Rydberg  atoms  and  Cooper  pair  boxes  [1 1 51. 1 1 611 .  The  molecular 
qubits  are  entangled  via  the  coupling  to  the  transmission  lines  rather  than  direct  dipole-dipole  in¬ 
teractions.  Again,  addressability  of  the  qubits  is  achieved  via  the  Stark  effect  by  means  of  local 
gating  of  an  electrostatic  field. 

Entanglement  is  a  major  ingredient  in  most  quantum  computation  algorithms.  It  is  among  the 


defining  features  of  quantum  mechanics,  with  no  classical  analog  I  17h19I1.  A  pure  state  of  a  pair 


of  quantum  systems  is  said  to  be  entangled  if  its  wavefunction  cannot  be  factored  into  a  product 
of  wavefunctions  of  the  individual  partners.  For  example,  the  singlet  state  of  two  spin-f  particles, 
-^=(|  tl)  -  I  IT))  is  entangled.  A  mixed  state  is  entangled  if  it  cannot  be  represented  as  a  mixture 
of  factorizable  pure  states.  The  allure  of  quantum  information  processing  has  recently  motivated 


studies  of  entanglement  for  a  variety  of  potential  qubit  systems  [5l  la  19, 


2014.3211 .  These  include 


one-dimensional  arrays  of  localized  spins,  coupled  through  exchange  interactions  and  subject  to 


an  external  magnetic  field  23]  and  analogous  treatments  of  trapped  electric  dipoles  coupled  by 
dipole-dipole  interactions  [321]  • 


However,  the  previous  studies  of  entanglement  of  electric  dipoles  have  not  adequately  consid¬ 
ered  how  the  external  electric  field,  integral  to  current  designs  for  quantum  computers  using  polar 
molecules,  affects  both  the  qubit  states  and  the  dipole-dipole  interaction.  For  the  simplest  case  of  a 
diatomic  molecule,  the  qubit  eigenstates  resulting  from  the  Stark  effect  are  linear  combinations 


2 


of  spherical  harmonics,  with  coefficients  that  depend  markedly  on  the  field  strength.  These  are  ap¬ 


propriately  termed  pendular  states  [33],  or  field-dressed  states  [34].  In  such  states,  the  orientation 


of  the  dipole  moment  has  a  broad  angular  range  (not  solely  along  or  opposed  to  the  field  direction 
as  are  spins  in  a  magnetic  field).  Likewise,  the  dipole-dipole  interaction  for  molecules  in  pendular 
states  is  much  different  than  that  for  dipoles  in  the  absence  of  an  external  field. 

Here  we  evaluate  entanglement,  as  measured  by  pairwise  concurrence,  for  the  prototype  case 
of  two  diatomic  polar  molecules  in  pendular  states,  ultracold  and  trapped  in  distinct  optical  lattice 
sites.  The  molecules  are  represented  as  identical  rigid  dipoles,  undergoing  angular  oscillations, 
a  fixed  distance  apart  and  subject  either  to  a  different  or  to  the  same  external  electric  field.  We 
examine  the  dependence  of  the  concurrence  on  three  dimensionless  variables.  The  first  governs 
the  energy  and  intrinsic  angular  shape  of  the  qubits  (when  the  dipole-dipole  interaction  is  switched 
off).  It  is  pE/B,  the  ratio  of  the  Stark  energy  (magnitude  of  permanent  dipole  moment  times 
electric  field  strength)  to  the  rotational  constant  (proportional  to  inverse  of  molecular  moment  of 
inertia).  The  second  variable  governs  the  magnitude  of  the  dipole-dipole  coupling.  It  is  Q/B.  with 
fl  =  (p2/r3),  the  square  of  the  permanent  dipole  moment  divided  by  the  cube  of  the  separation 
distance.  The  third  variable,  kBT /B,  is  the  ratio  of  thermal  energy  (Boltzmann  constant  times 
Kelvin  temperature)  to  the  rotational  constant. 

We  also  examine  an  aspect  related  to  but  distinct  from  entanglement.  The  operation  of  a  quan¬ 


tum  gate  [35]  such  as  CNOT  requires  that  manipulation  of  one  qubit  (target)  depends  on  the  state 


of  another  qubit  (control).  This  is  characterized  by  the  shift,  Aw,  in  the  frequency  for  transition 
between  the  target  qubit  states  when  the  control  qubit  state  is  changed.  The  shift  Aa>,  which  is  due 
to  the  dipole-dipole  interaction,  must  be  kept  smaller  than  the  differences  required  to  distinguish 


130  for 


among  addresses  of  qubit  sites.  Under  conditions  envisaged  in  the  proposed  designs  1H5], 
quantum  computing  with  trapped  polar  molecules,  Cl/B  <  10-4,  and  for  the  ground  eigenstate  both 
the  entanglement  and  frequency  shift  A  to  become  very  small.  For  CNOT  and  other  operations,  en¬ 
tanglement  needs  to  be  large,  but  can  be  induced  dynamically,  so  need  not  be  appreciable  in  the 
ground  eigenstate.  Yet  a  small  Aco  shift  can  only  suffice  if  the  resolution  is  high  enough  to  detect 
the  shift  unambiguously.  From  estimates  of  the  line  widths  of  transitions  between  the  pendular 
qubit  states,  we  find  it  an  open  question  whether  adequate  resolution  can  be  obtained  for  typical 
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candidate  diatomic  molecules. 


II.  ENTANGLEMENT  FOR  TWO  DIPOLES  IN  PENDULAR  STATES 


A.  Hamiltonian  terms  and  pendular  qubit  states 


The  Hamiltonian  for  a  single  trapped  linear  polar  molecule  in  an  external  electric  field  is 

p 2  2 

H  =  +  Vlrap(r)  +  BJ2  -  p  ■  £  (1) 

where  the  molecule,  with  mass  m,  rotational  constant  B  and  body-fixed  dipole  moment  p,  has 
translational  kinetic  energy  p2 /2m,  potential  energy  Vtrap  within  the  trapping  field,  and  rotational 
energy  BJ2  as  well  as  interaction  energy  p  £  with  the  external  field  8.  In  the  trapping  well,  at 
ultracold  temperatures,  the  translational  motion  of  the  molecule  is  quite  modest  and  very  nearly 
harmonic;  p2/lm  +  V,nip(r)  thus  is  nearly  constant  and  can  be  omitted  from  the  Hamiltonian.  There 
remains  the  rotational  kinetic  energy  and  Stark  interaction, 

Hs  =  BJ2  -  p£ cosd  (2) 


which  represent  a  spherical  pendulum  with  6  the  polar  angle  between  the  molecular  axis  and  the 
field  direction.  Figure  [H a)  displays  the  lowest  few  pendular  eigenenergies  [36]  for  a  'X  diatomic 
(or  linear)  molecule,  as  functions  of  p£/B.  These  are  labeled  with  the  familiar  quantum  numbers 


J,  M  that  specify  the  field-free  rotational  states.  However,  J  wears  a  tilde  to  indicate  it  is  no 
longer  a  good  quantum  number  since  the  Stark  interaction  mixes  the  rotational  states,  whereas  M 
(denoting  the  projection  of  the  J-vector  on  the  field  direction)  remains  good  as  long  as  azimuthal 
symmetry  about  8  is  maintained.  As  proposed  by  DeMille,  the  qubit  states  |0)  and  |1)  are  chosen 


as  the  lowest  M  =  0  pendular  states,  with  J  =  0  and  1,  respectively.  These  are  superpositions  of 
Yjo  spherical  harmonics, 


|0)  =  J^ajYj^cp),  |1)  =  J^,o  (G,(p)  (3) 

j  j 

Figure  [2]  plots  the  coefficients  as  functions  of  p£/B.  Figure  [3]  displays  the  angular  distributions 
of  the  pendular  qubit  states.  For  |0)  the  distribution  is  unimodal  and  as  p£/B  increases  the  dipole 
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orientation  increasingly  favors  the  direction  of  the  8-field  (at  6  =  0°).  For  |1)  the  distribution 
is  bimodal  because,  with  M  =  0,  the  dipole  is  rotating  perpendicular  to  the  J-vector,  which  is 
perpendicular  to  the  field  direction.  For  8  =  0,  the  dipole  orientation  is  equally  probable  in  the 
hemispheres  toward  {6  <  90°)  or  opposite  ( 6  >  90°)  to  the  field  direction.  As  pS/5  increases,  the 
pinwheeling  dipole  favors  the  opposite  hemisphere  because  there  its  motion  is  slowed  because  the 
Stark  interaction  becomes  unfavorable.  However,  when  pE/B  becomes  large  enough,  pinwheeling 
is  inhibited  and  converted  into  pendular  libration  about  the  field  direction,  so  the  dipole  orientation 
shifts  to  favor  the  toward  hemisphere. 

Adding  a  second  trapped  polar  molecule,  identical  to  the  first  but  distance  r12  apart,  introduces 
in  addition  to  its  pendular  term  the  dipole-dipole  coupling  interaction, 


Vd-d  = 


Bi  //2  —  3(//i  •  n)(//2  •  n) 
|ri-r2p 


(4) 


Here  n  denotes  a  unit  vector  along  r12.  In  the  presence  of  an  external  field,  it  becomes  appropriate 
to  express  Vd-d  in  terms  of  angles  related  to  the  field  direction.  As  shown  in  Appendix  A,  the  result 
after  averaging  over  azimuthal  angles  (that  for  M  =  0  states  are  uniformly  distributed)  reduces  to 


Vd-d  =  D(1  -  3cos2or)cos0icos@2 


(5) 


where  Q  =  p2/r3n,  the  angle  a  is  between  the  rl2  vector  and  the  field  direction  and  polar  angles  6\ 
and  92  are  between  the p\  and n2  dipoles  and  the  field  direction.  Until  later  (Sec.  IV),  we  consider 
the  external  field  magnitude  and  direction  to  be  the  same  at  the  sites  of  both  the  polar  molecules. 


B.  Entanglement  measured  by  pairwise  concurrence 


We  will  deal  with  the  entanglement  of  formation,  (£(p),  which  characterizes  the  amount  of 
entanglement  needed  in  order  to  prepare  a  state  described  by  a  density  matrix,  p.  (Henceforth,  we 


term  £(p) just  ’’entanglement”,  for  short.)  Wootters  1137. 


38]  has  shown  that  G(p)  for  a  general  state 


of  two  qubits  can  be  quantified  by  the  pairwise  concurrence ,  C(p),  which  ranges  between  zero  and 
unity.  The  relation  can  be  written  as 

«(p)  =  £(C(p))  (6) 
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where  £  is  given  by 


m  =  h 


1  +  Vl  -C- 


(7) 


with  h(x)  =  -x\og2x  -  (1  -  *)log2(l  -  x).  The  function  £(C)  increases  monotonically  between  zero 
and  unity  as  C  varies  from  0  to  unity.  The  concurrence  is  given  by 


C(p)  =  max  |0,  x[A\  -  -  xfh  ~  V^-} 


(8) 


where  the  T,’s  are  the  eigenvalues,  in  decreasing  order,  of  the  non-Hermitian  matrix  pp,  where  p  is 
the  density  matrix  of  the  spin-flipped  state,  defined  as 


P  =  (<Ty  <g>  Cry)p*(CTy  <8>  CTy)  (9) 

with  p*  the  complex  conjugate  of  p  and  cry  a  Pauli  matrix.  The  parent  density  matrix  p  is  taken 
in  the  basis  formed  by  combining  the  pendular  qubit  states;  for  a  pair  of  two-level  particles,  this 
comprises  the  four  state  vectors  ||00),  |01),  1 10),  1 1 1 )}. 

In  order  to  evaluate  thermal  entanglement,  we  need  a  temperature  dependent  density  matrix, 
p  =  exp (-/?//) /Z(T),  with  (3  =  1  /kBT  and  Z{T)  the  partition  function 

Z(T)  =  tr[exp(-/3//)]  =  ^  gie^E‘  (10) 

i 

with  Ef  the  ith  eigenvalue  and  g;  its  degeneracy.  Hence  the  density  matrix  can  be  written  as 

1  N 

(11) 

i 

where  |'Pi)  is  the  ith  eigenfunction.  From  the  density  matrix  p(T ),  we  can  obtain  the  reduced 
density  matrix  for  any  pair  of  dipoles  and  thence  evaluate  the  concurrence  at  any  temperature. 


III.  CONCURRENCE  OF  TWO  DIPOLES  IN  PENDULAR  STATES 

We  illustrate  the  calculation  of  pairwise  concurrence  for  =  2  dipoles.  The  Hamiltonian,  Hsl 
+  Hs 2  +  Vd-d,  when  set  up  in  a  basis  of  the  qubit  pendular  states,  ||00),  |01),  1 10),  1 1 1)},  takes  the 
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form, 


HS1+H 


S2 


W0  + W' 

0 

0 

0 

0 

Wo  +  w; 

0 

0 

0 

0 

Wi 

+  W' 

0 

,  0 

0 

0  Wj  +  w 

'  C0C' 

CoC'x 

CxC' 

CxQ 

-  3cos2a) 

CoC'x 

c0c; 

CXC[ 

CxC'0 

CiC0 

CiC'x 

,CxC'x 

CXC\ 

CiC'x 

Cic; 

(12) 


(13) 


where  W0  and  W\  are  the  eigenenergies  of  the  pendular  qubit  states  |0),  and  1 1 ),  in  the  absence  of 
the  dipole-dipole  interaction.  Primes  attached  to  quantities  for  the  second  dipole  indicate  that  the 
external  field  magnitude  may  differ  at  its  site  (although,  as  noted  above,  we  postpone  evaluating 
that  case  until  Sec.  IV).  In  Vd-d  the  basis  qubit  states  are  linked  by  matrix  elements  containing 
factors  arising  from  the  orientation  cosines  in  Eq.  ©;  these  are 


Co  =  (O|cos0|O);  Cx  =  <O|cos0|l>;  C,  =  <l|cos0|l> 


(14) 


C0  and  C\  are  the  expectation  values  of  cos 6  in  the  pendular  states  |0)  and  1 1),  respectively,  so 
represent  for  those  states  the  effective  dipole  moment  projections  displayed  in  FigOJb).  Cx  cor¬ 
responds  to  an  exchange  interaction  or  transition  dipole  moment  between  the  qubit  states.  Both 
the  Stark  eigenenergies  W,-  and  the  dipole-dipole  elements  Ck  are  functions  of  /dE/B.  As  seen  in 
FigEQb),  as  /dE/B  is  increased  C0  becomes  increasingly  positive,  whereas  C\  is  increasingly  neg¬ 
ative  until  about  /dE/B  =  2,  then  climbs  to  zero  at  about  /dE/B  =  4.9  and  thereafter  is  increasingly 
positive.  The  range  /dE/B  =  2  to  4  is  recommended  for  the  proposed  quantum  computer  designs 


rilL  Il6h:  within  that  range,  the  difference  in  the  effective  dipole  moments  of  the  qubits,  |C0  -  C i|, 
varies  only  modestly. 

If  the  dipole-dipole  interaction  is  omitted  (Q  =  0),  the  eigenvectors  of  Hsl  +  Hs2  are  simply 
Tj  =  |00),  ^2  =  2-1/2(|10)  -  |01)),  T'a  =  2~1/2(|10>  +  |01»,  W4  =  |11>,  corresponding  to  the 
eigenenergies  of  Eq.(12).  For  T/i  and  VP4,  which  are  obviously  nonentangled  states,  the  concur¬ 
rence  is  zero.  For  TS  and  'f'3,  which  exemplify  fully  entangled  states,  the  concurrence  is  unity; 
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these  are  termed  Bell  states  [170. 


When  the  dipole-dipole  coupling  is  included,  an  analytical  solution  to  obtain  eigenstates  is  only 
feasible  when  the  external  field  is  switched  off.  As  shown  in  Appendix  B,  in  that  limit  analytical 
results  can  be  obtained  for  each  step  in  evaluating  the  concurrence,  both  for  the  four  individual 
eigenstates  and  their  combination  in  the  thermal  concurrence.  As  seen  in  Fig.  1,  for  /uE/B  =  0,  the 
energy  terms  in  Eq.(12)  involve  merely  Wo  =  0  and  W\  =  2 B.  In  the  Vd-d  matrix  of  Eq.(13),  the 
cosine  matrix  elements  C0  and  C i  vanish  and  Cx  =  3~l/2;  thus,  the  only  nonzero  elements  occur 
along  the  antidiagonal  and  (for  a  =  90°)  are  just  QC|.  The  results  for  this  zero-field  limit  prove 
useful  in  interpreting  those  for  the  general  pendular  case. 

The  limits  with  D  =  0  and/or  /uE/B  =  0  motivate  setting  up  the  Hamiltonian  of  Eqs.(12)  and 
(13),  for  the  (unprimed)  case  with  the  same  external  field  at  both  dipole  sites,  using  a  basis  of  Bell 
states: 


111) +  |00)  111) -|00)  |10) +  |01)  |10) -|01) 


V2  ’  V2 

In  this  basis,  the  Hamiltonian  becomes 


V2  ’  V2 


Hsi  +  Hs2  — 


S2 


Vd_d  =  D(1  -  3cos 2 a) 


W+  W_  0  0 

W_  W+  0  0 

0  0  W+  0 

0  0  0  W+/ 

A+  B  C+  0 
B  A_  C_  0 
C+  C.  D+  0 


(15) 


(16) 


(17) 


v  o  o  o  b. , 

Where  W±  =  Wj  +  W,,  and  A±  =  ^(C2+C2)+C|,  B  =  ^(C2-C2),  C±  =  C^Q+C 0),  D±  =  C,C0±C\. 
This  makes  explicit  a  consequence  of  the  symmetry  between  the  (unprimed)  sites  [39].  In  the  Bell 
basis,  the  Hamiltonian  factors,  with  the  state  2  1/2(|  10)  - 101>)  in  a  1  x  1  block,  so  that  state  remains 
maximally  entangled  regardless  of  the  value  of  /jE/B  or  Q./B. 

Figure  H] plots,  for  [jE/B  =  0,  2  and  4.9,  the  eigenenergy  and  pairwise  concurrence  versus  Q/B 
=  0  to  6  for  the  four  eigenstates  of  the  two-dipole  system.  The  eigenstates  are  numbered  from 
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1  to  4  in  order  of  increasing  energy.  For  nE/B  =  0,  both  eigenstates  2  and  3  are  Bell  states, 
with  eigenenergies  Ej/B  =  2  -  (£2/65)  and  2  +  (£2/6 B),  respectively;  eigenstates  1  and  4  are 
also  entangled  (much  more  weakly)  by  the  dipole-dipole  interaction,  with  eigenenenergies  that 
shift  downwards  and  upwards  nonlinearly  with  increasing  £2/5,  respectively.  For  /uE/B  >  0,  the 
concurrences  increase  with  £2/5  for  eigenstates  1  and  4,  and  decrease  for  eigenstate  3.  By  virtue  of 
the  symmetry  imposed  factorization  noted  above,  eigenstate  2  retains  the  same  Bell  form  despite 
the  Stark  and  dipole-dipole  interactions  which  affect  its  energy,  and  its  concurrence  is  always  unity. 
For  small  Q/B  «  1,  eigenstate  3  also  becomes  independent  of  the  dipole-dipole  interaction 
and  coincides  with  eigenstate  2  in  both  energy  and  concurrence.  For  /rS /B  =  4.9,  as  seen  in 
Fig.  1(b),  the  C\  =  (l|cos0|l)  factor  that  appears  in  seven  of  the  matrix  elements  in  Eq.  (13) 
vanishes.  Consequently,  the  energy  of  eigenstate  4  then  becomes  independent  of  the  dipole-dipole 
interaction,  although  its  wavefunction  and  concurrence  do  not. 

Figure  [5]  shows,  for  [iE/B  =  0  and  2,  how  the  contributions  of  the  basis  states  to  each  of  the 
eigenstates  vary  with  the  strength  of  the  dipole-dipole  interaction.  This  illustrates  that  for  Q/B  « 
1  the  eigenstates  rapidly  approach  those  for  £2  =  0.  Indeed,  we  find  that  for  Q/B  <  0.04  the 
concurrences  for  eigenstates  1  and  4,  which  rapidly  become  the  same,  are  proportional  to  Q/B 
within  better  than  1%.  Thus, 

C\2-  K(x)[Q/B]  (18) 


where  the  proportionality  factor  K(x)  is  a  function  of  x  =  fiE/B.  At  the  zero-field  limit,  K (0)  =  1  /6. 
In  Appendix  C  we  describe  a  numerical  analysis  that  provided  an  accurate  approximate  formula, 


K{x)  =  Ay  + 


A2 


(19) 


1  +  exp[(x  -  x0)/Ax] 

This  is  plotted  in  Fig. 6  and  values  of  the  four  parameters  are  listed  in  Appendix  C. 

Figure  [7] displays  for  Q/B  =  0.1,  1  and  6  the  eigenenergies  and  concurrences  versus  ixE/B  from 
0  to  8  for  the  four  eigenstates.  As  the  dipole-dipole  interaction  increases  60-fold  over  this  range,  its 
effect  on  the  eigenstate  energies  is  relatively  modest,  whereas  the  concurrences  change  markedly, 
in  response  to  variations  in  eigenvector  compositions  such  as  illustrated  in  Fig.  0 

Figure  [8]  gives  a  contour  plot  of  the  thermal  pairwise  concurrence  derived  from  Eq.(ll)  as  a 
function  of  Q/B  and  kBT/B.  It  pertains  to  jiE/B  =  3;  we  found  that  normalizing  the  thermal 
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concurrence  to  its  value  for  T  =  0  and  Q/B  =  1  removed  most  of  the  variation  with  nE/B  from 
such  contour  plots.  For  T  =  0,  the  thermal  concurrence  coincides  with  that  for  the  ground  state, 
eigenstate  VP1.  However,  as  kBT/B  increases,  the  thermal  concurrence  decreases  and  is  always 
smaller  than  the  ground-state  concurrence.  This  may  seem  odd,  because  Eq.(l  1)  specifies  a  shift  in 
population  that  reduces  the  contribution  from  the  gound  state,  while  bringing  in  contributions  from 
the  excited  states.  The  eigenstates  2  and  3  then  populated  have  large  concurrence,  so  increasing 
temperature  might  be  expected  to  make  the  net  thermal  concurrence  become  larger  than  for  the 
ground-state,  rather  than  smaller.  The  source  of  this  behavior  is  indicated  by  the  analytic  solution 
obtained  in  Appendix  B  for  the  zero-field  limit, 

Ci2(T)  =  Cn(\)Px  -  C12(2)P2  -  C12(3)F3  -  C12(4)P4  (20) 


where  P,  =  ( l/Z)exp(-Ej/kBT )  with  Z(T )  =  X,  cxp(-E,/kBT ).  This  shows  that  the  excited  states 
indeed  reduce  the  thermal  concurrence,  an  effect  traceable  to  Eq.(8)  and  which  persists  even  for 
large  /uE/B. 

Another  striking  aspect  of  Fig.  8  is  that  the  concurrence  vanishes  along  and  outside  a  particular 
contour.  That  contour  defines  mutually  dependent  maximum  values  of  kBT /B  and  minimum  values 
of  Q/B  required  to  obtain  nonzero  concurrence.  When  Q / B  «  1,  we  find  that  a  modified  form  of 
Eq.(fT8l)  represents  the  thermal  concurrence, 


C12(T)  =  max{0,  K(x)[y  -  y0(x,z)]} 


(21) 


Here  x  =  nE/B\  y  =  Q/B:  and  z  =  kBT / B  is  the  scaled  temperature.  Fig.  [9]  gives  a  contour  plot 
of  y0  =  Qc/B,  the  critical  dipole-dipole  coupling  required  for  nonzero  concurrence.  Some  further 
details  are  included  in  Appendix  B. 

The  original  proposal  by  DeMille  and  kindred  papers  on  quantum  computing  with  trapped  polar 


molecules  [lj,  2, 


5|,llllOLll3 


16]  discuss  for  several  examples  the  range  of  experimental  conditions 


deemed  suitable  and  acceptable.  For  trap  temperatures  of  the  order  of  a  microkelvin  or  below, 
the  typical  values  of  kBT / B  are  a  few  times  1(T6,  so  indicate  that  only  ground-state  entanglement 
would  be  significant.  The  external  field  strengths  considered  are  typically  a  few  kV/cm.  The 
spacing  between  optical  lattice  sites,  r  =  A/2,  is  half  the  optical  lattice  wavelength.  The  optimal 
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choice  of  A  ranges  between 


to  0.3  microns,  depending  on  electronic  transition  frequencies  of 


the  molecules  to  be  trapped  [13].  From  these  parameters  and  molecular  data,  values  of  il/B  are 
small;  we  find  for  a  dozen  potential  candidate  molecules  values  ranging  between  4  x  10“6  (for 
KCs)  to  2  x  10~4  (for  Csl).  A  favorite  candidate  is  SrO  (ju  =  8.9 D,  B  =  0.33  cm-1,  A  =  1  micron), 
for  which  il/B  ~  10“5.  In  that  regime,  the  concurrence  is  simply  proportional  to  i'll B.  so  can  be 
easily  evaluated  from  Eq.(fT8l)  and/or  (I2TI)  without  use  of  the  rather  elaborate  prescription  outlined 
in  Eqs.(6-14). 


IV.  FREQUENCY  SHIFT  FOR  TWO  COUPLED  DIPOLES  IN  PENDULAR  STATES 


In  the  region  il/B  <  10  4,  the  concurrence  of  the  ground  eigenstate  is  very  small,  typically 


<  10  5.  However,  such  meager  entanglement  in  eigenstates  can  still  be 


computing,  as  demonstrated  with  NMR  versions  of  quantum  computers  [40].  The  key  aspect  is  that 


adequate  for  quantum 


although  entanglement  needs  to  be  large  for  some  quantum  computing  algorithms,  it  need  not  be 
appreciable  or  even  present  in  the  ground  eigenstate  of  the  system;  it  can  be  induced  dynamically 


during  operation  of  the  computer  [41].  Here,  for  the  polar  molecule  case,  we  consider  this  aspect. 


We  also  evaluate  an  eigenstate  property,  a  small  frequency  shift,  distinct  from  but  related  to  the 
pairwise  concurrence,  that  is  important  for  quantum  computing. 

The  need  for  selective  excitation  in  operation  of  quantum  logic  gates  [351,  421  43]  is  an  essential 
feature.  Taking  the  2-qubit  CNOT  gate  as  an  example,  its  operation  requires  that  manipulation  of 
one  qubit  (target)  is  perceptively  affected  by  the  state  of  the  other  qubit  (control).  In  our  case,  the 
qubits  are  pendular  states  that  can  be  accessed  by  microwave  transitions,  which  offer  high  spectral 
resolution.  As  resolution  has  a  crucial  role,  we  now  suppose  the  external  field  differs  enough  at 
the  two  dipole  sites  (denoted  unprimed  and  primed)  to  supply  distinct  addresses  for  the  sites  ( cf 
Fig.  1(a),  green  dashed  curve). 

Since  il/B  is  so  small,  we  first  omit  the  dipole-dipole  interaction  and,  as  illustrated  in  Fig.  [TO] 
consider  transitions  among  the  pendular  eigenstates  of  Eq.(12).  Although  in  this  limit  the  ground- 
state  concurrence  is  zero,  as  seen  in  Eq.(fT8l).  it  is  possible  to  generate  states  of  large  concurrence 


by  use  of  resonant  pulses  [41 


44].  Start  by  applying  a  pulse  resonant  with  the  transition  denoted 
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a>i,  between  |00>  and  |01),  which  has  energy  W[  -  W'y  Note  that  a>i  needs  to  be  well-resolved  from 
the  transition  a>3,  between  |00)  and  1 10),  which  has  energy  W\  -  W0.  The  separation  thus  comes 
from  the  different  values  of  the  external  field  at  the  two  sites  (plus  a  dipole-dipole  contribution, 
in  higher  order).  The  requisite  field  strength  difference,  8'  -  8,  can  be  readily  determined  from 
another  approximation  formula, 

W  -  Wo)/B  =Al+  A 2  (22) 

i  + (x/xo y 

by  comparing  results  for  x  =  /uE/B  and  x'  =  //£'/ B;  the  accurate  fit  obtained  (better  than  1%  except 
near  x  =  0)  is  displayed  in  Fig.  [[Hand  the  four  parameters  in  Eq.(l22l)  are  given  in  Appendix  C. 
The  amplitude  and  duration  of  the  (0\  pulse  can  be  adjusted  to  make  it  a  n/2  pulse,  which  will  put 
the  system  in  the  state  2_1/2(|00)  +  |01 )). 

Next,  to  complete  the  CNOT  gate,  apply  a  pulse  resonant  with  the  transition  co2  between  |01) 
and  |11).  This  needs  to  be  well-resolved  from  transition  a>3  between  |00)  and  1 10).  However,  in 
our  initial  approximation,  both  oj2  and  a>3  have  the  same  transition  energy,  Wj  -  W0.  Hence,  weak 
as  it  is,  the  dipole-dipole  interaction  is  seen  to  have  an  essential  role:  to  introduce  a  frequency 
shift,  A  a)  =  (jl>3  -  co2,  adequate  for  unambiguous  resolution.  If  that  is  fulfilled,  the  amplitude  and 
duration  of  the  co2  pulse  can  be  adjusted  to  make  it  a  n  pulse.  Thereby  the  system  will  be  put  in 
the  state  2-1/2(|00)  +  |1 1)).  This  result  of  a  CNOT  gate  is  to  first  approximation  a  Bell  state  (aside 
from  small  corrections  of  order  Q/B),  so  its  concurrence  will  be  near  unity.  It  is  not  an  eigenstate, 
so  will  evolve  with  time  but  in  principle  would  remain  nearly  fully  entangled  until  degraded  by 
other  interactions. 

If  now  the  dipole-dipole  terms  from  Eq.(13)  are  included  to  first  order,  we  obtain 


C0i  =  <01|#|01)  -  <00|#|00)  =  it;  -  W'0  +  n„C0(C;  -  Q  (23a) 

m2  =  <11|#|11)  -  <01|£|01>  =  W}  -  Wo  +  f \C;(Cj  -  Co)  (23b) 

co3  =  <10|#|10)  -  <00|tf|00)  =  Wi  -  Wo  +  QffC;(Ci  -  Co)  (23c) 

co4  =  <11|#|11)  -  <10|#|10)  =  w;  -  W'  +  f2ffCi(C;  -  C'0)  (23d) 

where  =  Q(1  -  3cos2cr).  Thus,  the  key  frequency  shift  is  given  by 

A  u)  =  co3  —  io2  —  (jj\  —  co\  =  f2ff(C[  —  c0)(c;  —  Cq)  (24) 
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For  given  Qa,  the  frequency  shift  Aco  depends  only  on  x  and  x',  which  determine  at  the  respective 
sites  the  difference  in  the  effective  dipole  moment  projections  Co  and  C\  along  the  external  electric 
field,  specified  in  Eq.(14).  To  provide  a  convenient  means  to  evaluate  Eqs.(23)  and  (24)  we  again 
fitted  our  numerical  results  to  obtain  accurate  approximation  formulas, 


Cq(x)  -  A  i  + 


i  +  (x/x0y 


(25) 


C i  (x)  -  Aq  + 


+ 


(26) 


1  +  exp[(x  -  xi)/Axi]  1  +  exp[-(x  -  x2)/Ax2] 

These  functions  are  plotted  in  Fig.  |TT]  together  with  C0  -  C\,  and  the  fitted  parameters  are  given 
in  Appendix  C. 

Since  for  small  Q./B,  both  the  concurrence  and  Aco  are  proportional  to  £l/B,  the  frequency  shift 
provides  an  equivalent  measure  of  entanglement.  When  the  8-fields  differ  at  the  two  sites,  Eq.(18) 
still  provides  a  very  accurate  approximation  for  Ci2(x,  x').  merely  by  replacing  the  proportionality 
factor  by  the  geometric  mean,  [K(x)K(x')\112 .  The  concurrence  (which  involves  Cx,  the  exchange 
interaction  term)  is  in  principle  different  from  Aco  but  both  have  about  the  same  magnitude.  The 
frequency  shift  is  much  more  relevant  for  quantum  computing,  because  Aco  is  directly  involved  in 
the  CNOT  gate. 

Also  important,  in  addition  to  the  pulse  shapes  which  affect  the  population  transfers,  are  the 
durations  of  the  resonant  pulses  required  to  resolve  coi  and  co2  from  coy,  these  must  satisfy  T31  >> 
l/|co3  -  coj  |  and  r32  »  l/|co3  -  co2 1.  For  r3i  the  lower  bound  usually  can  be  made  very  low, 
permitting  a  short  pulse  duration.  This  holds  because  a£  as  well  as  dipole-dipole  terms  contribute 
to  |  co3  -  ojf\,  which  thus  can  be  made  large  by  choice  of  the  8-field  gradient,  regardless  of  whether 
Qa  is  extremely  small.  In  contrast,  for  r32  the  separation  Aco  =  |co3  -  co2|  depends  only  on  the 
dipole-dipole  interaction.  The  smaller  Aco  is,  the  longer  the  co2  pulse  duration  has  to  be  in  order 
to  complete  the  CNOT  operation.  Although  larger  Aco  allows  a  shorter  pulse  duration,  Aco  must 
not  be  so  large  that  it  becomes  comparable  to  or  larger  than  the  addressing  shift  produced  by  a£, 
thereby  thwarting  correct  identification  of  the  qubits. 
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TABLE  I:  Illustrative  CNOT  Gate  Implementation13. 


/jE/B 

x  =  1 

x'  =  1.01 

x'  =  1.10 

x  =  3 

x'  -  3.03 

x'  =  3.30 

(Wi  -  W0)/B 

2.2709 

2.2759 

2.3218 

3.5614 

3.5831 

3.7789 

Co 

0.30165 

0.30404 

0.32487 

0.57922 

0.58149 

0.60051 

Ci 

-0.16467 

-0.16573 

-0.17461 

-0.16362 

-0.16150 

-0.14115 

o 

o 

1 

rs 

0.46632 

0.46977 

0.49948 

0.74284 

0.74298 

0.74165 

(o»i  -  m3)/fl 

4.99E-3 

5.09E-2 

2.17E-2 

2.17E-1 

A  lo/B 

2.19E-6 

2.33E-6 

5.52E-6 

5.51E-6 

Cn 

1.20E-6 

1.17E-6 

3.57E-7 

3.34E-7 

"For  D-„l B  =  10  5;  x  —  /yS /B,  x'  =  fj&'/B.  As  the  quantities  in  the  lowest  three  rows  are  functions  of 
both  x  and  x',  their  values  are  listed  in  the  x'  columns.  There  E-n  denotes  a  factor  of  10~". 


Table  I  provides  specific  numbers  pertaining  to  the  SrO  example.  From  Sec.III,  we  take  0.a/B  = 
10“5.  As  representative  8-field  values,  we  use  x  =  1  and  3  for  site  1  and  and  take  x'  higher  by  1% 
or  10%  for  site  2.  From  Eqs.(23),  the  transition  frequencies  a>\  =  a>4  and  <o2  =  m3  (in  units  of  B ) 
to  5  or  6  significant  figures.  The  frequency  difference  that  must  be  resolvable  for  the  first  step  of 
the  CNOT  operation,  tO\  -  m3,  is  approximately  just  A  A  IT  =  (Wj  -  IT(j)  -  (Wj  -  Wo).  From  Fig. 
mi  this  is  seen  to  grow  about  linearly  with  both  x  and  x'  -  x.  The  values  in  Table  I  (third  row  from 
bottom)  range  from  >  10“3  to  >  10_1  (in  units  of  B).  To  accommodate  more  dipole  qubits,  it  may 
be  desired  to  make  much  smaller  the  8-field  differences  between  sites;  steps  with  Ax  =  0.01% 
were  proposed  by  DeMille  [lj].  That  might  encounter  engineering  limitations,  but  in  principle 
the  proportionally  smaller  oj\  -  m3  difference  could  still  be  readily  resolved.  For  the  second  step 
of  the  CONT  operation  the  crucial  frequency  shift,  Am  =  a>i  -  m3  varies  only  modestly  with  x 
and  practically  not  at  all  with  x'  -  x.  The  values  of  Ato/B  in  Table  I  (second  row  from  bottom) 
range  between  2  and  6  x  1 0-6;  in  frequency  units,  this  range  is  20  to  60  kHz.  Smaller  still  are 
the  corresponding  values  of  the  concurrence  (bottom  row);  also  insensitive  to  x'  -  x  but,  in  accord 
with  Fig.  6,  varying  more  rapidly  with  x. 

Figure  d2]exhibits  for  both  the  A  to  shift  and  concurrence  the  variation  with  a,  the  angle  between 


14 


the  direction  of  the  electric  field  and  the  axis  between  the  dipoles.  This  dependence  enters  via  the 
factor  (1  -  3cos2ff)  in  the  dipole-dipole  interaction,  Eq.(5),  which  emerges  directly  in  the  A co  shift, 
Eq.(24),  and  by  a  more  complex  route  propagates  into  the  concurrence,  via  Eq.(13).  Tilting  the 
field  direction  to  make  a  =  54.73°,  the  ’’magic  angle”,  provides  a  simple  means  to  shut  off  the 
entanglement.  That  is  a  useful  option,  awkward  to  attain  in  other  ways  I7L 1 13M . 


V.  CONCLUSIONS  AND  PROSPECTS 


In  this  study,  our  chief  aim  has  been  to  examine  entanglement  of  polar  molecules  by  the  dipole- 
dipole  interaction  and  subject  to  an  external  electric  field,  for  the  prototype  case  of  two  diatomic 
or  linear  1  £  molecules.  This  required  use  of  qubits  that  are  pendular  states  comprised  of  sums  of 
spherical  harmonics.  We  focused  on  the  pairwise  concurrence  and  its  dependence  on  three  unitless 
reduced  variables,  involving  the  dipole  moments,  field  strength,  rotational  constant,  dipole-dipole 
coupling  and  temperature.  We  have  considered  a  wide  range  of  the  parameters,  to  map  general 
features  of  the  concurrence.  However,  for  conditions  envisioned  for  proposed  quantum  computers, 
the  dipole-dipole  coupling  is  weak  (Q/B  typically  of  order  10-4  to  1(T6)  and  the  concurrence 
becomes  very  small  (<  10  5).  For  that  weak  coupling  realm,  we  found  the  Aco  frequency  shift 
provides  an  equivalent  measure  of  entanglement,  directly  related  to  observable  properties  and 
hence  preferable  to  the  concurrence.  We  also  obtained  for  both  the  Aco  shift  and  concurrence  in 
the  weak  realm  simple  explicit  formulas  in  terms  of  the  reduced  variables. 

For  quantum  computing  a  crucial  issue  is  whether  Aco  is  large  enough  to  enable  the  co3  transition 
to  be  reliably  distinguished  from  co3  (and,  equivalently,  col  from  co4).  For  typical  candidate  polar 
molecules,  this  requires  resolving  transitions  separated  by  only  tens  of  kHz.  That  would  not  be 
feasible  in  conventional  molecular  spectroscopy.  Under  ordinary  gas  phase  conditions,  transitions 
between  molecular  rotational  or  pendular  states  have  line  widths  of  the  order  of  a  few  100  kHz 


[45].  For  ultracold  molecules  trapped  in  an  optical  lattice,  line  widths  may  be  much  narrower. 


Collisional  broadening  is  eliminated  and  at  microkelvin  temperatures  Doppler  broadening  is  also 
quenched  (as  trap  conditions  are  in  the  Famb-Dicke  regime).  It  is  encouraging  that  for  ultracold 
atoms  extremely  narrow  line  widths  have  been  attained  by  exploiting  ’’magic”  optical  trapping 
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conditions  that  are  expected  to  be  at  least  in  part  applicable  to  molecules  [46] .  At  present,  however, 


no  data  have  been  reported  on  line  widths  for  rotational  transitions  of  ultracold  molecules  trapped 
in  an  optical  lattice  and  subject  to  an  external  electric  field.  In  view  of  the  small  size  of  Aw,  it 
is  important  to  obtain  such  data  to  assess  the  resolution  attainable,  since  motion  within  the  traps, 
coupling  to  lattice  fields,  and  inhomogeneity  of  the  external  field  may  introduce  appreciable  line 
broadening. 

We  have  sought  to  glean  pertinent  evidence  from  electric  resonance  spectroscopy  of  molecular 
beams,  as  the  beams  are  collision  free  and  transitions  are  observed  in  an  external  electric  field 


(”Rabi  C-field”)-  For  BaO,  both  A M  =  0,  J  =  0  — »  1  transitions  in  the  microwave  region  [47]  and 


A  J  =  0,  \M\  —  0  — >  1  transitions  in  the  radiofrequency  region  [48]  have  been  observed,  in  fields 


ranging  from  ~  200  -  500  V/cm.  For  the  radiofrequency  transitions,  line  widths  were  only  about 
2  kHz,  consistent  with  just  the  dwell  time  in  the  C-field.  But  for  the  microwave  transitions  the 
widths  are  much  larger,  45  kHz;  this  is  attributed  both  to  the  higher  frequency  of  the  transitions 
and  to  experimental  conditions  that  render  more  significant  Doppler  broadening  and  nonunifor¬ 


mity  of  the  field,  especially  in  the  entrance  and  exit  fringe  regions  [49].  The  Doppler  and  dwell 


time  contributions  are  not  relevant  to  inferring  what  might  be  expected  for  trapped  BaO  (or  SrO). 
Broadening  by  inhomogeneity  of  the  external  field  is  relevant  but  depends  very  much  on  experi¬ 
mental  particulars.  The  transitions  of  interest,  depicted  in  Fig.  [Hi  occur  in  the  microwave  region 
and  involve  Stark  fields  typically  ten-fold  larger  than  used  in  the  electric  resonance  spectroscopy, 
so  the  line  widths  might  be  significantly  broadened  due  to  field  inhomogeneity.  These  observations 
do  not  permit  firm  conclusions  about  the  resolution  issue,  but  it  decidedly  poses  an  experimental 
challenge. 

This  discussion  pertains  only  to  the  choice  of  qubits  we  have  considered,  pendular  states  of 
linear  polar  molecules,  which  involve  transitions  that  change  J  but  not  M.  The  resolution  issue 
motivates  examining  other  choices  for  qubits.  For  instance,  states  with  the  same  J  but  different 
M  could  be  used.  Other  options,  particularly  use  of  hyperfine  or  nuclear  spin  states  instead  of 
pendular  states,  have  been  suggested  as  means  to  reduce  sources  of  decoherence  II lL l2L  1 1 3l.  1 1 611 .  As 
yet,  the  size  of  A co  for  any  qubit  choice  other  than  that  used  in  this  paper  remains  to  be  determined. 

We  intend  to  extend  the  treatment  developed  here  to  other  choices  for  qubit  basis  states  as 
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well  as  to  larger  numbers  N  >  2  of  dipoles.  In  preliminary  work  on  linear  and  planar  arrays  of 
dipoles  up  to  N  =  8,  we  find,  as  expected,  the  maximum  pairwise  concurrence  occurs  for  next- 
neighbor  dipoles,  although  that  for  non-nearest  ones  is  significant.  Also  in  prospect  ij^an  analogous 


treatment  of  the  proposed  coupling  of  polar  molecules  via  microwave  strip-lines  [16].  There  the 


entangling  interaction  differs  from  direct  dipole-dipole  interactions,  but  again  the  proposed  qubits 
are  pendular  states. 

Another  pendular  variant  inviting  attention  is  use  of  polar  symmetric  top  molecules.  The  |0) 
and  |1)  qubits  can  be  selected  as  |  J,K,M)  =  |1, 1,  -1)  and  1 1, 1,  +1),  which  are  degenerate  in  the 


field-free  limit  and  thus  have  a  first-order  Stark  effect  [45]].  Even  in  a  weak  electric  field,  these 


states  are  strongly  oriented  along  and  opposed  to  the  field,  with  equal  and  opposite  projections. 
Moreover,  the  effective  dipole  moments  do  not  depend  on  the  field  strength  so  low  fields  can  be 
used  if  necessary  to  reduce  line  broadening,  without  the  penalty  imposed  by  quenching  of  effective 
dipoles  that  would  occur  for  the  second-order  Stark  effect.  At  first  blush,  the  symmetric  top  option 
appears  to  be  disallowed  because  transitions  between  the  M  =  -1  and  +1  Stark  states  violate  the 
selection  rules,  AM  =  0  or  A M  =  ±1.  But  the  prohibition  is  not  absolute.  Because  the  optical 
lattice  perturbs  cylindrical  symmetry  about  the  field,  M  is  not  strictly  a  ’’good”  quantum  number, 
so  the  AM  selection  rule  is  relaxed.  Moreover,  if  the  molecule  contains  an  atom  with  nuclear  spin 
I  >  1/2,  and  hence  an  electric  quadrupole  moment,  transitions  with  A M  =  +2  become  allowed. 


For  instance,  a  deuterium  nucleius  (/  =  1)  makes  AM  =  +2  transitions  facile  in  Stark  spectra  [50]. 


Other  symmetric  top  options  for  qubits  are  inversion  doublets  (e.g.  in  NH3)  or  internal  rotation 
states  associated  with  hindered  torsional  motion  (e.g.  CH3CF3);  these  offer  strong  dipole-allowed 
transitions. 

previous  studies  of  entanglement,  both  for  polar  molecules  [7, 


03 


27 


32]  and  for  magnetic  spins 


310,  have  considered  primarily  domains  where  the  concurrence  is  large  (>  0.1),  and  have 


focused  on  means  to  tune  the  entanglement  to  attain  such  domains.  For  polar  molecules,  that 
requires  Q /B  >  1.  Recently,  it  was  suggested  that  such  large  Q /B  could  be  attained  for  dipole 


arrays  by  exploiting  nanotraps  with  lattice  spacing  of  the  order  of  only  10  nm  [32].  However, 


as  emphasized  in  Sec.  IV,  for  quantum  computing  large  entanglement  in  the  ground  eigenstate 
is  not  required.  Indeed,  reducing  the  array  spacing  so  markedly  would  strongly  foster  inelastic, 
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spontaneous  Raman  scattering  of  lattice  photons  and  hence  induce  decoherence  [lj,  2, 


ij.  Such 


considerations  make  small  rather  than  large  Q./B,  and  consequently  weak  rather  than  strong  en 


tanglement  in  the  ground  eigenstate,  actually  preferable  for  quantum  computing  [37],  provided 
resolution  of  the  A  to  shift  can  be  attained. 


ACKNOWLEDGEMENTS 

We  are  grateful  to  the  Army  Research  Office  for  support  of  this  work  at  Purdue  and  to  the 
National  Science  Foundation  (CFIE-0809651),  the  Office  of  Naval  Research  and  Hackerman  Ad¬ 
vanced  Research  program  for  support  at  Texas  A&M  University.  We  thank  David  DeMille  and 
Seth  Lloyd  for  instructive  discussions  clarifying  the  role  of  entanglement  in  quantum  computing; 
Christopher  Ticknor  for  calling  attention  to  kindred  aspects  of  ’’dressing”  dipole-dipole  collisions 
in  an  external  electric  field;  and  William  Klemperer  for  elucidating  field  effects  in  electric  reso¬ 
nance  spectroscopy  and  how  to  evade  the  usual  Stark  transition  selection  rules. 


APPENDIX  A:  DIPOLE-DIPOLE  INTERACTION 


The  angular  dependence  of  the  dipole-dipole  interaction,  given  in  Eq.@,  is  usually  expressed 


as 


<D ij  =  cos/3  -  3cosy;Cosy;- 


(Al) 


where  j3  is  the  angle  between  dipoles  //,  and  Hj\  angles  y,  and  y;  specify  the  orientation  of  the 
dipoles  with  respect  to  the  vector  r,;  between  them.  Ordinarily,  it  is  natural  (and  done  in  all 
textbooks)  to  express  cos/3  in  terms  of  the  angles  y  together  with  the  azimuthal  angles  <f>r  about  the 
Tjj  axis.  Thus,  use 


cos/3  =  cosy,cosy;-  +  siny;siny7cos(0n  -  <prj) 


(A2) 


which  when  combined  with  the  -3cosy,cosy;  term  gives  the  familiar  expression  [51]].  In  the  pres¬ 


ence  of  the  external  electric  field,  we  need  to  recast  <f>(/  in  terms  of  angles  0,  and  6j  that  specify  the 
orientation  of  the  dipoles  with  respect  to  the  direction  of  the  external  electric  field.  Therefore,  we 
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use 


cos ft  =  cosdjcosdj  +  sindjsind jCos(cf>Ei  -  <f>Ej )  (A3) 

cosy  =  cosdjcosa  +  sin0,sinQ'COs(0£/  -  (pEr)  (A4) 

where  the  4>e  are  azimuthal  angles  about  the  field  vector  E  and  a  is  the  angle  between  the  field 
vector  and  the  rtj  vector.  The  azimuthal  factors  can  be  expressed  as 

cos (<pEi  ~  <f>Ej)  =  cos(pEiCOS(/)Ej  +  &w0Eisinfej  (A5) 

As  M  =  0  states,  which  do  not  depend  on  the  </>  angles,  are  chosen  as  the  qubit  basis  states,  in 
evaluating  matrix  elements  of  <D;j  between  these  states  the  integrations  over  d(pid(pj  (from  0  to  2/r) 
eliminate  all  terms  involving  the  </>£  angles.  The  net  result  is  simply 

=  (1  -  3cos2a)cosd,cos0j  (A6) 

The  effect  of  integrating  over  the  6  angles  is  just  to  replace  in  Vdd  the  dipole  moments  q,  and  jid  by 
their  effective  values,  fi(cosB).  The  effective  dipole-dipole  interaction  hence  reduces  to 

Vd-d  =  £2(1  -  3cos2a)(cos0,)(cos0y)  (A7) 

with  Q  =  HiUj/r?.  as  a  convenient  scale  factor. 


APPENDIX  B:  ZERO-FIELD  CASE 


For  ju£ /B  =  0,  the  Hamiltonian  matrix  reduces  to  diagonal  terms  from  Eq.(12)  and  antidiagonal 
elements  from  Eq.(13),  and  the  pendular  qubit  basis  states  become  simply  |0)  =  Too  and  |1)  =  Ti0. 
The  form  of  the  Hamiltonian  makes  it  equivalent  to  that  for  the  Ising  model  for  a  system  with  two 


qubits  [52].  Diagonalization  of  the  Hamiltonian  yields  the  eigenenergies  and  eigenvectors  given 


in  Table  II  as  explicit  functions  of  £2/5. 
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TABLE  II:  Zero-Field  limit  for  N  =  2  dipoles. 


i 

Eigenenergies,  E,-/B 

Wavefunction,  T, 

Cl2 

i 

2-2(1  +£2)1/2 

— r==?(|H}  -  a+|00» 

y\+a+ 

2  «+ 

I+q-2 

2 

2(1-0 

^(|10>  -  |01» 

1 

3 

2(1+0 

^(|10>  +  |01» 

1 

4 

2  +  2(1  +  £2)1/2 

— p=(|1  1)  -  a_|00» 

"\/l+0C 

2|o-| 

l+oi 

where  a±  = 


i±(i+<r2) 


1/2 


/£,  with  £  =  Q./6B 


The  density  matrix  for  eigenstate  1,  the  ground-state,  is 


p(  1)  =  I'PtX'Pt 


1  +  a; 


That  for  eigenstate  2  is 


P(  2)  =  P2)<^2I  =  +- 


a:  0  0  — ctj 


0  0  0  0 

0  0  0  0 

-QT+  0  0  1 


^00  0  oN 


0  1-10 
0-110 
0  0  0  0 


(Bl) 


(B2) 


The  p(3)  matrix  differs  from  p(2)  by  having  +1  in  place  of  each  -1;  the  p(4)  matrix  differs  from 
p(l)  by  having  a_  in  place  of  a+.  As  these  density  matrices  pertain  to  only  two  dipoles,  they  need 
not  be  reduced  further. 

Obtaining  the  density  matrices,  p(z),  for  the  spin-flipped  states,  defined  in  Eq.  ©,  involves 
shuffling  the  rows  and  columns  of  p(z)  in  accord  with  |00)  ++  |1 1)  and  |01)  ++  1 10).  This  gives 


P(l)  = 


1  +  on 


1  0  0  -a. 


0  0  0  0 
0  0  0  0 


-a+  0  0a 


(B3) 


+ 
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and  thus  the  product  matrix  is 


p(  1)P(1)  = 


(i  +  4) 


2  a2  0  0  -a\  ^ 


0  0  0  0 
0  0  0  0 


-2a+  0  0  2  a2  j 


(B4) 


2 

The  eigenvalues  of  this  matrix  are  A\  =  4<r;/  (l  +  cr+)  ,  T2  =  T3  =  T4  =  0.  From  Eq.  ([8]),  the 
concurrence  is,  Ci2(l)  =  2a+/(l  +q7  ).  Similarly,  we  find  the  concurrences  for  the  other  eigenstates, 
given  in  Table  I. 

To  evaluate  the  thermal  concurrence,  we  need  to  set  up  the  thermal  density  matrix 


where 


with 


Then  we  find 


p(T)  =  YjzxV(-EilkBT)\'¥i)(x¥i\ 


1=1 


^00/ 
0  b  d  0 
0  d  b  0 
g  0  0  c 


a- 


1  +  ai 


-P\  + 


a 


1  +  ac 


b  =  ~(P2  +  Ps) 

1  1 
P 1  +  - 7P4 


1  +  a 


1  +  az 


d  =  -(P3  -  P2) 


O', 


«  = 


1  +  a: 


-P 1 


cr_ 


1  +  ar_ 


= 


exp  (~Ej/kBT) 


p(T)p(T)  = 


ac  +  g-  0  0 

0  b2  +  d2  2  bd 

0  2  bd  b 2  +  d2 


lag 

0 

0 


leg 


0 


0  ac  +  g2 


(B5) 


(B6) 

(B7) 

(B8) 

(B9) 

(B10) 

(Bll) 


(B 1 2) 
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and  obtain  the  eigenvalues, 


d^lV^-g)2,  A2  =  (b-d)2,  A3  =  (b  +  d)2,  d4  =  (V^  +  g)2  (B 13) 

Hence  from  Eq.  d8]),  we  obtain  the  thermal  concurrence, 

Cn(T)  =  max  {0,  -2  (b  +  g)}  (B 14) 

With  the  Cl2(i)  of  Table  I,  this  gives  Eq.(20)  of  the  text, 

CniT)  =  C12(1)P,  -  C12(2)P2  -  C12(3)P3  -  C12(4)P4  (B15) 

When  Q/B  «  1,  the  ground  state  concurrence  becomes  Ci2(l)  — »  D/6 B,  in  accord  with 
Eq.(fT8l)  of  the  text,  whereas  Ci2(4)  — >  D/6 B  and  C]2(2)  =  C]2(3)  =  1.  Provided  that  also 
kBT /B  «  1,  a  first-order  expansion  of  Eq.(B15)  gives 

C\2(T)  «  C12(l)  -  e2-  63,  (B16) 

where  e2  =  P2/P\  «  1;  e3  =  P3/P2  «  1  and  P4/P1  «<  1.  Then 

Cn(T)  *  K{ 0)  [Q/B  -  Q.JB] ,  (B 17) 

which  has  the  form  of  Eq.(21)  of  the  text  with  K{ 0)  =  1/6  and 

y0(T)  =  Q.JB  =  (e2  +  e3)/K(  0)  =  l2exp(-2B/kBT)cosh(Q/3kBT)  (B18) 

This  result  for  the  zero-field  case,  although  not  useful  in  practice,  illustrates  how  the  excited 
states  are  involved  in  creating  a  temperature  dependent  minimum  level  of  dipole-dipole  coupling, 
QJB,  that  is  required  to  have  nonzero  thermal  concurrence. 

Figure  IT31  shows  a  contour  plot  of  C]2(T)  for  the  zero-field  case,  derived  from  Eq.(B14).  It  is 
qualitatively  quite  similar  to  Fig.  8  for  the  pendular  case,  over  the  same  wide  range  of  kBT / B  and 
a/B. 

APPENDIX  C:  REDUCED  VARIABLE  FORMULAS 

In  order  to  find  a  proper  reduced  variable  formula,  three  steps  are  needed:  (1)  calculate  enough 
sample  points  to  define  well  the  exact  curve;  (2)  find  a  function  with  adjustable  parameters  that 
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enables  fitting  those  points;  (3)  evaluate  the  parameters  using  a  non-linear  regression  method. 
For  our  curve  fitting  we  use  the  Levenberg-Marquardt  Algorithm  [53],  also  called  ”Chi-square 
minimization”.  Chi-square  is  defined  as: 


x1  =  Yj[Yi-m,e)f  (Cl) 

i=  1 

where  xt  and  V,  are  the  independent  and  dependent  variables  for  the  ith  (i  =  1 ,  2,...,n)  sample  points 
of  the  exact  curve;  0  are  the  parameters  to  be  fitted.  The  Levenberg-Marquardt  algorithm  iteratively 
adjusts  the  parameters  to  get  the  minimum  chi-square  value,  which  corresponds  to  the  best  fit.  The 
input  data  for  fitting  Eqs.(19, 22, 25, 26)  comprised  our  numerical  results  for  the  pendular  case, 
over  the  ranges  x  =  /uE/B  =  0  to  8.  Tables  III  -  VI  list  the  optimal  values  found  for  the  parameters 
and  95%  percent  confidence  intervals.  At  x  =  0,  the  Eq.(19)  fit  gives  K(Q)  =  0.17103,  different 
slightly  from  the  exact  zero-field  limit,  K(0)  =  1/6.  Likewise,  at  x  =  0  the  Eq.(25)  and  (26)  fits 
give  Co  =  0.005  and  C i  =  0.00072  rather  than  the  exact  value  of  zero.  The  critical  point  for  C\  to 
change  sign  is  x  =  4.902,  whereas  Eq.(26)  gives  C\  =  -0.00025,  slightly  different  from  zero. 


TABLE  III:  Values  of  the  parameters  for  Eq.(fl9l). 


Parameters 

Pendular 

Field-free 

CP' 

Ai 

0.01092 

0.00221 

±0.0003 

a2 

0.21953 

0.24779 

±0.006 

Xo 

0.96578 

0.74035 

±0.05 

Ax 

0.97429 

0.86072 

±0.03 

a  95%  confidence  interval;  values  listed  are  maximum  found  for  the  2  curves 
shown  in  Fig[6]  Both  R2  values  are  around  0.9981. Similarly  accurate  results 
are  found  when  Ea.(fl8l)  is  generalized  for  different  E-fields  at  the  dipole  sites 
by  replacing  K(x )  by  [K{x)K{x')]111. 
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TABLE  IV:  Values  of  the  parameters  for  Eq.(l22l ). 


Parameters 

Values 

CP 

Ai 

12.42379 

±0.0533 

A2 

-10.47646 

±0.0534 

Xo 

8.77516 

±0.0534 

p 

1.5867 

±0.00527 

a  95%  confidence  interval.  R 2  =  0.9999 

TABLE  V:  Values  of  the  parameters  for  Eq.(l25l). 

Parameters 

Values 

CP 

Ai 

0.84855 

±0.00145 

a2 

-0.84355 

±0.00180 

x0 

1.6339 

±0.00508 

P 

1.2459 

±0.00539 

a  95%  confidence  interval.  R2  =  0.99994 

TABLE  VI:  Values  of  the  parameters  for  Eq.(l26l). 

Parameters 

Values 

CP 

Ao 

-0.75212 

±0.0323 

Ai 

1.04192 

±0.0336 

a2 

1.14092 

±0.0325 

xi 

-0.16241 

±0.0224 

X2 

3.1232 

±0.124 

AX\ 

0.90544 

±0.0136 

AX2 

2.76286 

±0.0496 

“95%  confidence  interval.  R2  -  1 
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For  convenience,  we  give  formulas  for  the  three  unitless  ratios,  evaluated  with  customary  units: 
l iE/B  =  0.0168  p(Debye)£(kV/  cm)  /  B(cm~l); 

Q/B  =  5.04  x  10-9 ir(Debye)/r3 (microns) 
kBT/B  =  0.695  T(K)/B(cm~l). 
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FIG.  1:  (Color  online)  Stark  states  for  a  polar  diatomic  molecule  in  a  electronic  state  [36D,  as  functions 


of  qE/fi,  with  p  the  permanent  dipole  moment,  8  the  field  strength,  B  the  rotational  constant,  (a)  Eigenen- 
ergies,  W,  and  (b)  Matrix  elements  of  orientation  cosines;  see  Eq.(14).  States  used  as  qubits  (red  curves) 
are  labeled  |0)  and  |1).  In  the  field-free  limit,  |0)  correlates  with  the  ./  =  0,  Mj  —  0  and  |1)  with  the  /  =  1, 
Mj  =  0  rotational  states.  Dashed  curve  (green)  in  (a)  shows  energy  for  transition  between  qubit  states, 
A W  =  W\  -  Wq',  that  in  (b)  shows  Co  -  C\,  difference  between  effective  dipole  moments,  projections  of  the 
molecular  dipole  on  the  field  direction  for  pendular  states  |0)  and  |1). 
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FIG.  2:  (Color  online)  Coefficients  of  sums  of  spherical  harmonics  for  pendular  states  |0)  and  |1),  see  Eq.(3). 
Dashed  curve  for  1 1 )  indicates  the  coefficient  of  0  is  negative. 
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FIG.  3:  (Color  online)  Angular  distributions  of  the  |0)  and  |1)  pendular  states  for  values  of  p£//f  between  0 
and  8. 
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FIG.  4:  (Color  online)  Eigenenergies,  numbered  1,  2,  3,  4  in  order  of  increasing  energy,  and  pairwise 
concurrences  of  the  eigenstates  for  two  dipoles  as  a  function  of  the  dipole-dipole  coupling  constant  D./B  for 
three  values  of  the  reduced  electric  field  strength,  /uE/B  =  0,  2  and  4.9. 
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FIG.  5:  (Color  online)  Eigenvectors  of  the  four  eigenstates  for  two  dipoles  as  a  function  of  the  dipole-dipole 
coupling  constant  D./B  =  0  to  6  for  p£//>’  =  0  (dashed  curves)  and  2  (solid  curves). 
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FIG.  6:  (Color  online)  The  K(x)  function  of  Eq.(18);  solid  curves  show  exact  result  (red)  and  fitted  function 
(blue)  of  Eq.(19)  that  pertains  to  pendular  qubit  basis,  see  Eq.(3).  For  comparison,  dashed  curve  pertains  to 
field-free  basis  with  |0)  =  Fq,o>  U)  =  Ti,o  (Cf  Table  III,  Appendix  C). 
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FIG.  7 :  (Color  online)  Eigenenergies  and  concurrences  for  the  four  eigenstates  for  two  dipoles  as  a  function 
of  reduced  variables,  /rS /B  for  electric  field  and  Q./B  for  dipole-dipole  coupling. 


34 


FIG.  8:  (Color  online)  Contour  plot  of  thermal  pairwise  concurrence  for  two  dipoles,  for  /rS /B  =  3.  For 
0  <  Q/B  <  1,  the  maximum  concurrence  Cnimax)  =  0.0473,  occurs  at  T  =  0,  Q./B  =  1.  The  plot  displays 
normalized  contours.  Within  each  colored  band,  the  variation  of  Cn/Cxiimax)  is  0.1;  thus  the  normalized 
concurrence  in  the  right  most  band  (red)  ranges  from  0.9  to  1;  in  the  next  band  (orange  red),  from  0.8  to  0.9, 
etc.  A  striking  feature  is  the  large  region  (uncolored)  where  Cn  =  0.  There,  entanglement  does  not  occur 
unless  the  dipole-dipole  coupling  exceeds  a  critical  value  dependent  on  the  temperature. 
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FIG.  9:  (Color  online)  Contour  plot  displaying  yo(x,z)  term  in  Eq.(21)  vs.  x  =  /uE/B.  Within  each  colored 
band,  the  range  of  >>o  is  0.001;  thus  in  the  lowest  colored  band  (magenta).  )>o  ranges  between  0  and  0.001; 
in  the  highest  colored  band  (red),  \>o  is  between  0.009  and  0.01. 
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FIG.  10:  (Color  online)  Schematic  energy  levels  for  qubit  pendular  eigenstates  of  N  =  2  dipoles,  in  absence 
of  dipole-dipole  interaction,  thus  corresponding  to  Eq.(12).  Qubit  basis  states  shown  at  left,  eigenenergies 
at  right.  Pairs  of  transitions  involved  in  CNOT  operation  are  indicated:  a>i  transfers  dipole  2  from  |0)  to 
|1)  with  dipole  1  remaining  in  |0);  then  a>2  transfers  dipole  1  from  |0)  to  1 1 )  with  dipole  2  remaining  in 
|1).  Analogously,  the  same  result  could  be  reached  by  m3  followed  by  to  4.  Transition  energies  (including 
dipole-dipole  terms  to  first-order)  are  given  in  Eqs.(23). 
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FIG.  11:  (Color  online)  Comparison  of  exact  results  (blue  curves)  with  fitted  approximation  functions 
(dashed  red  curves)  for  properties  governing  transitions  among  qubit  states,  Eqs.(23):  pendular  energy 
difference,  (W\  -  Wq)/B ,  cosine  expectation  values,  Co  and  C\  and  their  difference,  Cq-  C\  \  cf  Eqs.  (l22l 

HU  |26]>. 
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FIG.  12:  (Color  online)  Frequency  shift  A co/B  (left  panels)  and  ground  state  concurrence  Cn  (right  panels) 
as  functions  of  a ,  the  orientation  angle  of  the  electric  field.  Curves  are  shown  for  Q./B  =  0.2  to  1.0  with 
luE/B  =  1  or  3. 
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FIG.  13:  (Color  online)  Contour  plot  of  thermal  pairwise  concurrence  for  field-free  case,  prepared  in  same 
format  for  comparison  with  Fig.  8  for  the  pendular  case.  Here,  Cu(max)  =  0.1644  at  T  =  0,  D./B  =  1. 
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